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ABSTRACT. 

We probe to what extent not only the counts of the faint blue galaxies and their 
redshift distribution but also their ^-resolved luminosity functions can be explained in 
terms of binary merging (aggregations) . 

We present a dynamical theory of such interactions. On this basis we find that 
"minimal aggregations" taking place within large scale structures and triggering star- 
bursts yield rates and timing such as to explain the observations of the local, flat lumi- 
nosity function and of its progressive rising and steepening for redshifts out to z ~ 1. 
Correspondingly, we predict faint blue counts still rising up to ~ 28 — 29, redshift 
distributions shifting toward larger z with increasing rag; in addition, we predict an 
upturn of the faint end of the luminosity function more pronounced in clusters than in 
the field. 

We propose that our picture provides the differential dynamics missing in the canon- 
ical hierarchical clustering theories. This reconciles with the observations the steep 
luminosity functions predicted at high z by such theories. 



Subject headings: cosmology - galaxies: evolution - galaxies: formation. 
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1. INTRODUCTION 



The observations of faint blue galaxies (FBGs) provided statistically significant counts 
in works by Koo (1986), Tyson (1988), and Cowie et al. (1988). Then an excess became 
apparent, namely, by a factor 5 at magnitudes B ^ 24 and ~ 7 at fainter ones, over 
the values expected from normal local galaxies filling a critical Friedmann-Robertson- 
Walker universe. Over the last decade the growing data base concerning FBGs kept 
challenging the interpretations. 

Various contributions to that excess have been discussed, the simplest being based 
on improved measurements of the local luminosity function (hereafter local LF). These 
include the correction of the incompleteness by a factor 1/2 in the previous samples 
of local galaxies, as pointed out by Lilly et al. (1995) and by Ellis et al. (1996), 
and the consideration of the upturn at Mb > —16 from a flatter Schechter shape 
pinpointed by Marzke, Huchra, & Geller (1994). Both these improvements in the local 
baseline increase, albeit moderately, the expected counts, since the latter are given 
by the convolution along the line of sight of the LFs weighted with the cosmological 
volumes. 

Meanwhile, it became apparent that toward faint magnitudes blue galaxies make 
up a fraction of all galaxies which increases both with decreasing luminosity and with 
increasing redshift z (Koo & Kron 1992; Colless et al. 1993). Thus, some sort of galaxy 
evolution must be at work. 

Passive luminosity evolution driven by the shift with age of stellar populations 
would provide the observed blue counts only if helped by the larger volume contributed 
by widely open cosmologies with Q D <^ 1 (see Pozzetti, Bruzual, & Zamorani 1996). A 
similar result also holds in the red band down to K = 24 (Cowie et al. 1994; Djorgovski 
et al. 1996), though the excess is considerably smaller here. Such cosmologies, however, 
are out of tune with other evidence (summarized, e.g., by Ostriker & Steinhardt 1995) 
favoring a flat universe of higher density, i.e., fi ^ 0.3 with a cosmological constant 
A D = 1 — fi . 

On the other hand, a constraint was added with the advent of large spectroscopic 
samples (Cowie, Songaila, & Hu 1991; Lilly 1993; Colless et al. 1993; Glazebrook 
et al. 1995) that provided redshift distributions N(z) with a median value as low as 
(z) ~ 0.5-0.6 for the identified objects with B ^ 24. The low value of the redshift where 
N{z) is peaked rules out strong luminosity evolution in a flat universe; rather, for 
z ^ 0.5, the distribution resembles the behavior provided by passive-evolution models 
with sophisticated color mixes (Koo 1990). The latter, on the other hand, would imply 
a very substantial high-z tail unless very specific recipes are implemented (see Gronwall 
& Koo 1995), which imply mismatches at lower z (see Cowie et al. 1996). 

However, recently Lilly et al. (1995) and Ellis et al. (1996) have provided direct 
evidence of considerable evolution in the resolved luminosity functions over the range 
z ~ 0.1 — 1. This concerns only the blue objects, while the red ones appear to be already 
"in place" by z ~ 1, and quite similar to typical L* galaxies of the local kind mainly 
undergoing passive change. The observed evolution appears to be strongly differential, 
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being faster for the lower luminosities measured. 

Two main interpretations have been proposed for the onset of such a blue, later 
lineage of galaxies. One (Babul & Rees 1992; Babul & Ferguson 1996) holds these 
to be very small and numerous lumps in which gas collapse and star formation have 
been delayed by the high level of the UV background prevailing before z « 1; shortly 
thereafter, the photoionizing UV flux might decline sharply enough for this constraint 
to be lifted. Following the initial short starburst, these objects would rapidly fade away. 

The other interpretation stems from the widespread local evidence of starbursts 
associated with galaxy-galaxy interactions (among the vast literature see the recent 
papers by Zabludoff et al. 1996; Hutchings 1996 and references therein); this maintains 
that many more, if smaller, blue starbursts are triggered at increasing depths by tidal 
interactions taking place mainly with or between numerous dwarfs (Lacey & Silk 1991). 
Some such interactions are likely to end up in galaxy-galaxy aggregations, i.e., binary 
merging, also causing considerable number evolution. A phenomenological treatment 
of these effects has been outlined by Guiderdoni & Rocca-Volmerange (1991) and by 
Broadhurst, Ellis, & Glazebrook (1992). 

A dynamical approach to the kinetics of aggregations has been given by Cavaliere & 
Menci (1993), motivated to complement the standard hierarchical clustering theory by 
its intrinsic inadequacy in two specific respects: it predicts too many small objects, and it 
envisages too little dynamics to erase them; so by itself it would set and preserve a steep 
luminosity function. A reconciliation between excess counts and a locally flat LF may 
be sought in terms of star formation processes (see White 1994; Frenk, Baugh, & Cole 
1995), but it still falls short of the full goal even on tuning several parameters. Cavaliere 
& Menci (1993) instead, on the basis of the dynamics of aggregations, predicted that 
N(M, z) and hence N(L, z) should flatten out with z decreasing in the range z ^ 1. 

We are encouraged to pursue our dynamical approach further not only by the 
recent data concerning the ^-resolved luminosity functions (which we now compare 
with) but also by some very recent lines of additional evidence. These include the 
following: a sizable fraction of FBGs at z ~ 1 show dynamical signs of ongoing or 
recent interactions, or the presence of blue satellites (Koo et al. 1996; Van den Bergh 
et al. 1996); some FBGs look like spirals with a red bulge (Koo et al. 1996) and 
peripheral starbursts; widespread post-starburst galaxies mark the action of pervasive 
galaxy-galaxy interactions (Zabludoff et al. 1996); measurements of clustered redshifts 
have shown the presence of at least some galaxy concentrations (Le Fevre et al. 1994; 
Koo et al. 1996). 

Here we will start from first principles and build up a dynamical theory of in- 
teractions containing no free parameters once one specifies the environment where the 
interactions take place. Our main point will be that many such interactions can occur 
even at moderate redshifts (z ~ 1) within large scale structures that offer the best com- 
bination of volume and density contrast. Such interactions, with their mass-dependent 
cross section £(M), naturally yield a differential evolution. This affects the dwarf galaxy 
numbers yielding considerable density evolution, while the associated starbursts affect 
all blue luminosities yielding a milder luminosity evolution. In terms of observables, 
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we find counts still rising at very faint magnitudes, redshift distributions peaked at a 
value shifting from z ~ 0.5 at B ^ 24 toward z ~ 1 for fainter magnitudes, and resolved 
luminosity distributions which that flatten out from z ~ 2 toward a local Schechter-like 
shape. 

The plan of the paper is as follows. In §2 we present our dynamical framework 
to compute the mass functions evolving by aggregations. In §3 we discuss the strong 
dependence of such evolution on the galaxy environment. We compare three kinds of 
environments: the homogeneous "field" ; the bound, virialized structures; and the large- 
scale structures, which in fact will dominate the observables. In §4 we describe how from 
the evolutionary mass distribution we compute the observables: the resolved luminosity 
distribution N(L,z), the blue counts jV(mg) and the redshift distribution N(z). In §5 
we give our results and predictions for these observables. In the final section we give 
and discuss our conclusions. 



2. THE MASS FUNCTION FROM BINARY AGGREGATIONS 



2.1 The dynamics 

A dynamical theory for galaxy evolution directly concerns the mass function N(M,t). 
The governing equation will be of the general form 

8N N (2.1) 



dt t ' 

where r _1 is an average evolutionary rate. 

In the simplest case when the evolution of a galaxy is not affected by surrounding 
companions, the rate will not depend on the mass function itself. But this produces 
self-similar N(M, t) with faint-end slope unchanged in time. A well known instance is 
constituted by the Press & Schechter (1974) formula, corresponding to the choice of r 
discussed by Cavaliere, Colafrancesco, & Scaramella (1991) and by Kitayama & Suto 
(1996). 

Hence, we consider the case when the evolutionary rate r -1 does depend on the 
number of surrounding masses that interact and eventually aggregate. For a total 
density n g (t) of galaxies, this corresponds to considering r _1 oc n g (SI/), where E is 
the cross section for aggregation and V is the relative velocity of the two aggregating 
masses; the average runs over the galaxy velocity distribution. The requirement of total 
mass conservation (see Lucchin 1988) requires a positive and a negative term on the 
right hand side of equation (2.1), and easily leads to the following Smoluchowski form 
for the kinetic equation: 



dN 1 



dt 2 



I dm K(m\ m — m', t) N(m', t) N(m — m', t) 
Jo 

/•OO 

-N(m,t) / dm' K( m,m ,t)N(m' ,t) . (2-2) 
Jo 
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Here the kernel K(m,m',t) = n g (t)(EV) corresponds to the interaction rate for each 
pair of masses m, m', and contains the physics driving the evolution of the distribution 
N(m,t). The latter is written in its comoving form; the masses m = M/M* are nor- 
malized in terms of M*, the mass corresponding locally to the standard characteristic 
luminosity L* (see Peebles 1993). 

Each galaxy aggregation comprises two stages: first, the dark matter halos coalesce; 
second, the baryonic cores may eventually coalesce. The former process requires grazing, 
weakly hyperbolic encounters for which the cross section is given by (Saslaw 1986, 
Cavaliere, Colafrancesco, & Menci 1992, hereafter CCM92) 



E = e (^-) 7r ( r - + r -') 



I G (m + m' 



V 2 



with r m = r* m 1 ^ 3 , v m = v* m 1 ^ 3 



(2.3) 

where and r* are the velocity dispersion and dark-halo radius of the present M* 
galaxy. We use M* = 5 10 11 h~ l M© 1 ; r* is conservatively taken to be 40 h~ l kpc (see 
Persic, Salucci, & Stel 1996), and v* = 270 km/s (three dimensional). The function 
£(V/vm) describes the decreasing efficiency of the aggregations for increasing relative 
velocities. Following the results from N-body simulations (see Richstone & Malumuth 
1983) we will take e = 1 for V < Sv m , and just e = for larger values; we shall discuss 
the effect of such a cutoff in §5 and §6. The second term (proportional to G) in the 
cross section describes the focusing effect of gravity, an important addition to the purely 
geometrical cross section proportional to n only for large masses and in the range 
V < V3v m that prevails in long-lived poor groups (see CCM92). 

The second stage is constituted by the possible coalescence of the baryonic cores. 
Sufficient conditions for this to occur have been discussed extensively in terms of dy- 
namical friction of one core against the background provided by the merged halos (see, 
e.g., Kauffmann, White, & Guiderdoni 1993). The associated scale tf is at most a factor 
of 5 longer than the halo crossing time when halo aggregations are effective, namely, for 
grazing encounters with V < 3v m , as discussed above. In all such cases, we find that tf 
does not exceed the halo aggregation time r which drives the evolution of N(M, t) after 
equation (2.1). Note that such a condition is a sufficient and conservative one, as dis- 
cussed by Kauffmann et al. (1993), and as confirmed by detailed simulations (including 
hydrodynamics) of single collisional events and of the ensuing aggregation process (see 
Barnes & Hernquist 1991, Barnes 1992). 

The solutions of equation (2.2) depend on the setting of the initial conditions, not 
only as a starting point of the integration procedure, but also - in keeping with the 
quadratic nature of the right-hand side - for establishing the amplitude of the kernel. 
We discuss these points in turn. 

2.2 Initial conditions 

The initial conditions are set up at the time ti n corresponding to a redshift Zi n . Because 



1 We parameterize the Hubble constant using H Q = 100 h km/s/Mpc. 
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we think of aggregations as complementary to the standard hierarchical clustering (for 
which see Peebles 1993), we focus on initial mass functions having the Press & Schechter 
(1974) form 

N(m) dm = 2 "_^ c 0(m) dm , (2.4) 

where 4>(m) = m~ 2+a exp(— b 2 5 2 m 2a /2). Here 5 C ~ 1.68, b~ l ~ 1 fixes the amplitude, 
and a = (n + 3)/6 is related to the spectral index n ~ — 2 of the power spectrum for the 
density fluctuations from which galaxies form (see White 1994). 

In the same spirit, the typical galactic mass at the initial time M*(ti n ) is identified 
with the characteristic mass M c (ti n ) provided by the standard hierarchical theory; when 
normalized to its present value [e.g., M c (t ) = 0.6 x 10 15 hr 1 M & for = 1], this reads 
(Lupini 1995) 

M c (t in ) 



M c (t a ) 



l/3a 

(2.5) 



for a density parameter fi D and for a cosmological constant satisfying the canonical rela- 
tion A G = 1 — o characteristic of most inflationary scenarios (see Ostriker & Steinhardt 
1995), with a = tt H q t /2fll /2 . 

In the following, we will consider mainly 6=1, and n = const ~ —2.5 on galactic 
scales. 

The initial mass function (2.4) not only provides the starting condition for equation 
(2.2), but also affects the numerical value of the kernel at the initial time t in . In fact, 
we obtain 

n g (Un) = Cn ° Pco ~^= 'WJf 1 ) J ^ m ^ dm ' ( 2 - 6 ) 

integrated from M*/500 (the smallest mass considered); here p co is the critical cosmo- 
logical density at the present epoch, and C makes explicit the contrast factor at t in of 
the environment where interactions take place. 

Note that assuming the Press & Schechter expression (2.4) and the above values 
for the parameters is not mandatory; because of the very negative values of the spectral 
index n at galactic scales, the exponent of the power-law in <p(m) is always close to -2, a 
generic feature of the mass distributions also found in other variants of the hierarchical 
clustering scenario (see, e.g., White 1994). On the other hand, shortly after ti n the 
integro-differential nature of equation (2.2) causes the solutions to lose memory of the 
detailed shape of the initial conditions. 



2.3 Interaction rates 

The numerical value of the interaction kernel K = n g (t){T, V) determines the timescale 
for the aggregation process. Computing n g as described above, and substituting fiducial 
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values for r* and for the initial relative velocity V^ n , we obtain an expression for K that 
we conveniently write in the form 

K(m, to', t) = 3 10~ 4 H 1 (m 2/3 + m' 2/s ) 

The function f(Cl , A ,t) = n s (t) V(t)/n ff (ti n ) V(tj n ) describes the change of the kernel 
after the initial ti n (depending, we recall, on the choice of A ), and includes the 
evolution of the relative velocities and that of the density p(t) within the environment. 
The factor t~ x makes explicit the unit of time we use throughout the computation. 

The initial density n g (t in ) is given by equation (2.6). Its value depends on the 
cosmo logical scenario under consideration; e.g., for an initial time £j n corresponding to 
a redshift Zi n = 3 and for a CDM spectrum (n ~ —2.2 at galactic scales) with = 1 
and C = 1, we obtain M*{t in ) = 5 x 10 11 h~ l M , so that n g (t in ) 1.8 h 3 (1 + z in ) 3 w 
1.2 x 10 2 /i 3 Mpc" 3 down to M*/500. 

The solutions of equation (2.2) are unique (McLeod 1962) once we have specified 
the initial conditions (through eqs. [2.4] — [2.6]) and the aggregation rate (equation 
[2.7]). We stress that these in turn depend on the following quantities: C, U n , Vi n , 
f(Q ,X ,t); their values describe the environments where aggregations are operating, 
as will be discussed next. 



l + ^(m 2 / 3 +m' 2 / 3 ) 



x 



3. ENVIRONMENTS 



The kinetic equation (2.2) is naturally written in terms of t. In the following we will 
derive in terms of z the quantities entering this equation and the resulting observables; 
the conversion from t to z for flat universes may be found, e.g., in Peebles (1993). 

3.1 The homogeneous "field" 

If any such description is realistic - which we doubt in view of the impressive evidence 
of large scale structures recalled below - the contrast factor should be taken as C = 1; 
the density n g oc (1 + z) 3 decreases after ti n in keeping with the general expansion; the 
relative velocities, given an initial value, scale as V oc (1 + z). Thus, the aggregation 
rate changes with z according to 

M^) 4 ' (3 ' 3) 

where zi n is the redshift corresponding to the initial time £j n . 

The corresponding rapid decrease of /(£) with time makes the aggregation rate 
rapidly negligible, so that the evolution of the mass function after equation (2.2) will 
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become mild at most. Nevertheless, for comparison purposes we compute it using a 
value for Zi n typical of the formation era of galaxies, according to equation (2.5) which 
in turn depends on fi Q , A Q and on the perturbation spectrum. For a scale-free spectrum 
with n = —2.2 we obtain Zi n = 3.5 when = 1 and A D = 0; and Zi n = 4.5 when 
fi = 0.2 and A G = 0.8. The initial Vi n = V a (1 + Zi n ) is given in terms of the local value 
V for the pairwise galaxy velocity dispersion. 

The latter was measured first by Davis & Peebles (1983) (see also Tormen et al. 
1993) who found V a ~ 340 km/s. Recently Marzke et al. 1995 combined the data from 
the CfA Redshift Survey and the Southern Redshift Survey and obtained a larger value 
540 km/s; however, when the Abell clusters with R > 1 are removed from the sample, 
the pairwise "velocity dispersion of the remaining galaxies drops to V Q = 298 km/s". 
We adopt the value 300 km/s. Actually, such a value for V Q is conceivably due to the 
present perturbed gravitational potential; so our scaling back as 1 + z in this subsection 
will optimistically maximize the effect, which will be shown to be inadequate anyway. 



3.2 Virialized structures 



Following the standard hierarchical clustering, virialized structures like groups or clus- 
ters as environments for galaxy interactions may be characterized simply in terms of 
the redshift z at virialization. Then, the density contrast is ~ 180, the galaxy density 
scales as n g (z)/n g (0) = (1 + z) s , and the virial velocity dispersion as 

71-1 

V(z)/V = (l + z)^+^ (3.4) 



for Q, = 1; here the spectral index takes on values in the range n = —1.5, —1.3 appro- 
priate for the scale of groups and clusters. Thus, the average aggregation rate would 
scale as the simple product of these two ^-dependencies. However, we have to con- 
sider the limited survival time of the structure before it is reshuffled into larger and 
less dense condensations. Such time is evaluated to be r p = 3t/2 on average (Cav- 
aliere, Colafrancesco, & Scaramella 1991; Lacey & Cole 1993), which corresponds to 
t p = St in (1 + z in f/ 2 /2 (1 + zfl 2 for ft = 1 and A D = 0. Including the suppression 
factor due to limited survival, the effective aggregation rate is r _1 r p /t D , and in full 
reads 



2rt + l 
n + 3 



n g (E V) r p /t = n g (E V)(z in ) | ^ ±±±- ) . (3.5) 

2 t \ L + Zi 



This is equivalent to considering the scaling 



2 n. + 1 



^ - \TTt) ' • (3 ' 6) 

and an effective contrast 

C = 270 t in / 1 . (3.7) 
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We start our computation at a redshift typical for formation of groups of galaxies, 
Zi n = 2.2 for an index n = —1.3 at the cluster scale. With this value, the initial Vi n 
is evaluated from equation (3.4) on the basis of a (three-dimensional) value V Q = 1270 
km/s for clusters of richness 1, to obtain Vi n = 480 km/s. The corresponding effective 
contrast from equation (3.7) is C « 45 for = 1. 



The large scale structures (LSSs; among the vast literature see De Lapparent, Geller 
& Huchra 1992; Sahni & Coles 1995; Vettolani et al. 1996) constitute an interesting 
environment for interactions to take place in; over the field they have the advantage of a 
larger density, and over groups and clusters the advantage of lower velocity dispersions. 
However, they are quantitatively less understood as yet from both the observational 
and the theoretical points of view. 

In the following, we shall use the semi-analytical approach proposed by Do- 
roshkevich et al. (1995) based on the Zel'dovich approximation and supported by nu- 
merical N-body simulations (see again Doroshkevich et al. 1995, and references therein). 
In such a picture, caustics form first in the matter field at Zf ~ 4 (for = 1), then 
LSSs build up by accretion of matter falling onto them. In this stage aggregations are 
not effective, because most galaxies have either closely parallel motions or large relative 
velocities. 

Aggregation can be effective in a subsequent stage when the thickness of the pan- 
cakes approaches half of their typical separation. This is the epoch Zi n when we start 
our calculations. We evaluate it taking up the results by Doroshkevich et al. (1995), for 
= 1. For z ^ 3 the mean comoving distance C between LSSs increases for decreasing 
z after 



where y = (1 + z)/(l + Zf) and r c depends on the power spectrum. The LSS thickness 
in this stage turns out to be 



According to the discussion above, we identify z in with the epoch when i\/C\ is 1/2. 
This yields z in « 2.3, independent of r c . 

The final stage for z < Zi n constitutes a comparatively quiet era during which 
the behavior of the galaxy relative velocities flattens out (Menci & Valdarnini 1993) to 
the observed value Vi n ~ 300 km/s (see Davis & Peebles 1983; Tormen et al. 1993). 

During this era, the scale of voids evolves as £2 ~ [\A + [(1 + zt) /(l + z)] 2 — l] ^ 2 
(Doroshkevich et al. 1995), with z t ~ 0.2. From such an expression we find the comoving 
thickness of the structures to scale approximately as (1 + z) 1 / 2 so that the contrast 



3.3 Large scale structures 




-1 



(3.8) 




(3.9) 
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behaves like C(z) oc (1 + z)^ D ~ 3 ^ 2 , where D is the dimensionality of the LSS (D = 2 
for sheet-like structures, and D = 1 for filaments). In terms of the function f(z) defined 
in equation (2.7) this translates to 

r -i (3+D)/2 

f(z)=[(l + z)/(l + z in )\ . (3.10) 

If the local contrast is ~ 5, consistent with observations (Geller & Huchra 1989; 
Ramella, private communication) and with numerical simulations (Brainerd and Vil- 
lumnsen 1991; Menci & Valdarnini 1993), then from the scaling of C(z) we infer an 
initial contrast C = 5(1 + z in )( D ~ 3 ^ 2 , with the value C 3 for D = 2. 



In this section we have defined the quantities which enter equation (2.2) for the 
computation of the evolving mass function in different environments. A summary of 
the resulting numerical values is given in Table 1. Next we describe how, from the mass 
function, we compute the observables like the luminosity function, the galaxy number 
counts and redshift distributions. 



4. FROM N(M, t) TO OBSERVABLES 



Our basic assumption in deriving the observables is that interactions also trigger 
star formation (see, e.g., Zabludoff et al. 1996, Hutchings 1996, and references therein). 
Specifically, we conservatively assume that in the range of parameters where aggrega- 
tions are effective, they also induce and sustain star formation at a rate s(z) equal to 
the rate of dynamical interactions r~ 1 (z) acting on the available gas (see Lacey & Silk 
1991; Broadhurst, Ellis, & Glazebrook 1992). 

The associated luminosity for the dwarf galaxies is given by the relation L oc 
Mg as /T, considering that r exceeds the dynamical time of a single galaxy and con- 
stitutes the limiting factor. Moreover, the deposition of energy from Supernovae will 
be initially balanced by cooling, but then for shallow potential wells it will go predom- 
inantly into gas expulsion. In this regime, even more directly than in White & Frenk 
(1991), the relation M gas ocMw^ is obtained. Since r _1 has the same M and t (or z) 
dependencies as given for K by equation (2.7), we obtain 

L/L* = (M/M*y and L*(z) oc f(z,X ,Q ) , (4.1) 

where r\ = 4/3. Such a value applies when the cross section (2.3) is purely geometrical 
and is consistent with the observational results by Kormendy (1990). From equation 
(4.1) the luminosity function N(L) = N[M(L)} dM/dL is obtained. 

Since our aim is to isolate the effects of dynamical evolution, we deliberately simplify 
colors and fc-corrections to the bones. So we identify the prompt luminosity L with that 
in the blue band (see Broadhurst et al. 1992) , and in the blue spectral region at moderate 
z we assume a neutral fc-correction (see Lilly, Cowie, & Gardner 1991 and Koo & Kron 
1992). We take an absolute characteristic magnitude M B * = —20 (see Lilly et al. 1995; 
Ellis et al. 1996). 



10 



The counts and the redshift distributions are found by convolving the comoving 
N(L, z) with the cosmological volume V, to obtain 



dz — (l + zf N[L(m b , z)] -j— dm B , (4.2) 

where the relation between L and ms is given, e.g., by Guiderdoni & Rocca-Volmerange 
(1991), and the cosmological volume for a flat universe with o + A D = 1 is given, e.g., 
by Peebles (1993). The range of integration over z extends out to the maximum z from 
which the magnitude tub can be observed. 



5. RESULTS 



In this section we present our results for the mass distribution, the z-depen- 
dent luminosity function N(L, z), the counts iV(ms) and the redshift distributions N(z) 
in different magnitude ranges. These observables are derived from galaxy evolution 
driven by aggregations acting in the different environments described in §3. 

The plots we present are obtained through the following steps: 

i) For a given galaxy environment we compute the aggregation rate given by equation 
(2.7); the parameters C, ti n (or Zj n ), Vi n and the function f(0, o ,X o ,t) are discussed in 
§3, and summarized in Table 1. 

ii) We numerically integrate equation (2.2) from the initial time ti n to the present time 
t a . The time step is set at At = (t — t o )/500 while the mass step is M*/500. 

iii) Following the equations recalled in §4 we compute N(L, z), iV(ms) and N(z). 
The results we obtain are compared with the observations concerning the faint blue 
counts, the ^-distributions, and the ^-resolved LFs. 



5.1 The homogeneous "field 



For comparison, we put on record the corresponding results, although this is not ex- 
pected to be an environment conducive to many late binary aggregations due to the 
rapid decrease of density and velocities. 

Figure la shows that such is indeed the case for O = 1 (A D = 0). Because of the fast 
decrease of f(Q ,X ,t) with time (see Table 1), the interaction rate rapidly becomes 
small, so that N(M, t) hardly evolves from the initial shape. The corresponding LFs 
are shown in Figure lb. Any small evolution that occurs takes place at high redshifts, 
so that no appreciable change can be observed in the range z ^ 1. The faint counts in 
Figure lc are well fitted, which is not not surprising considering the steep LFs at all 
redshifts. But the ^-distributions (Figure Id) are peaked at z ~ 1 even at relatively 
bright magnitudes. 

Thus, aggregations of galaxies uniformly distributed in a FRW universe are un- 
interesting, and clearly fail to account for the observed evolution in the LFs out to 
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One might expect that the plateau (actually the inflection) characteristic of the 
cosmic scale factor in a Lemaitre universe (A 7^ 0) could provide an era of roughly con- 
stant density when the binary aggregation processes could proceed for a while. However, 
in a flat universe such processes are suppressed after the inflection by the accelerated 
expansion. The net effect is still close to the case O = 1, as can be seen from Figures 
2a — 2d computed for fi G = 0.2 and A D = 0.8. 

5.2 Virialized structures 

A more favorable environment for galaxy aggregations is constituted by the bound, 
virialized structures like groups and clusters of galaxies. Here the large density contrast 
C (though convoluted with the limited survival time discussed in §3.2), and the milder 
time decrease of the interaction kernel /(0 o , A D , t) yield a sustained evolution of N(M, t) 
(see Figure 3a) compared to that found for the "field" . In Figures 3b — 3d we show the 
LFs, the counts, and the ^-distributions expected if all galaxies resided within virialized 
structures for a considerable fraction of their lifetime. It is found that the evolution of 
N(L,z) consists of a rapid flattening in the sub-L* range, with relatively little change 
in the range from z ~ 1 to the present. 

It must be noted that in clusters, in spite of the larger value of the contrast C, the 
aggregation rate is suppressed not only by the limited survival time (see §3.2), but also 
by the increase with decreasing z of the velocity dispersion entering the cutoff e(V/v m ) 
in equation (2.3). The latter circumstance has the effect of decreasing the gravitational 
term in E, and even of shrinking to zero the whole cross section of the smaller galaxies 
with the lower internal velocities; so at the faint end the mass distribution will hardly 
be changed from the initial steep shape. The net result will be a LF with an upturn 
toward the faint end (see Figure 3e). The balance between the effect of the sharp, 
mass-dependent decrease of the cross section and the general speeding up driven by the 
large effective contrast is a delicate one, and deserves an aimed study motivated by the 
LFs recently found (see Driver et al. 1994) to be appreciably steeper in local clusters 
than in the "field" . 

5.3 Large scale structures 

A result in a way intermediate between the two extremes above is provided by the 
environment constituted by the LSSs, where observations show most galaxies to reside 
not only locally (Geller & Huchra 1989), but also at increasingly larger redshifts (see 
Broadhurst et al. 1990; Schectman et al. 1995; Vettolani et al. 1995; Carlberg et al. 
1996). For galaxies inside sheet-like structures (our case D = 2) the parameters entering 
the aggregation rate are summarized in Table 1, and yield the evolution of N(M,t) 
shown in Figure 4a. The corresponding counts and redshift distributions (Figures 4c, 
4d) agree with the observations. In particular, they agree with the spectroscopic redshift 
distribution by Glazebrook at al. (1995) if L oc M 4 / 3 holds, and yield an enhanced high- 
z tail for a stronger dependence of L on M, e.g., L oc M 15 . The luminosity functions 
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(Figure 4b) show a clear flattening from zrjI toward a Schechter-like, local shape with 
—1.1 in agreement with the data. 

Note that an upturn at faint luminosities should also be expected here; the mecha- 
nism is similar to that discussed at the end of §5.2, but the effect is considerably milder 
(and shifted toward fainter luminosities) because the relative velocities inside LSSs are 
and remain smaller than those in clusters. 

Defining a characteristic luminosity by the ratio of the second to the first moment 
of N(L, z), we find this evolves by about 1 from z = 1 to z = 0, as also illustrated in 
Figure 4a. Actually, the intrinsic luminosity change L*(z) oc (1 + z)( 3+D ^ 2 (see eqs. 
[4.1] and [3.10]) is nearly balanced by the change following approximately (1 + z) -1 " 5 , 
due to characteristic mass after the relation L oc M 4 / 3 . The net result is a moderate 
change (about 1 Mag) of all luminosities, competing with the large change in number for 
setting the shape of N(L, z). The corresponding M/L ratio evolves (for an galaxy) 
as M^z)- 1 / 3 (1 + z )- (3+D V 2 which yields M*/L* ~ (1 + z)-( 2+D V 2 . 

6. CONCLUSIONS AND DISCUSSION 

We have derived for FBGs a luminosity function that changes considerably for 
increasing z as a result of number evolution and a milder luminosity evolution. In 
particular, in the observed range out to z ~ 1 the luminosity function N(L, z) rises and 
steepens from a local logarithmic slope about —1.1 to a value close to —1.6 at redshift 
w 0.8, as observed (see Figure 4b, and below for a discussion). 

This is derived from a dynamical theory of the galaxy mass function based on mass- 
conserving binary aggregations. The resulting evolution of N(M, t) is far from being 
self-similar; indeed, it is strongly differential with mass, causing for smaller galaxies 
stronger number changes and shift of the mass scale. This constitutes an important 
complement to the canonical hierarchical clustering, which in fact contains too little 
dynamics to account for the observed evolution of the luminosity function, as recalled 
in the Introduction. 

The actual strength of the evolution induced by aggregations depends on the envi- 
ronment where these take place. In CCM92 we stressed that in long-lived groups strong, 
gravitationally focused interactions of large galaxies occur, leading to the formation of 
cD-like mergers. Here, we find that in the large-scale structures — which comprise most 
local galaxies, are being observed out to z ~ 1, and in a critical universe are expected 
to loom out at redshifts z ~ 2 — small-small and small-large aggregations are favored. 
Their strength is such as to match - with no special effort at optimization - the existing 
data for the resolved N{L, z) and for the integrated but deeper observables like counts 
and magnitude distributions (see Figures 4a to Ad). 

A simple way of summarizing the interaction (and starburst) rates derived in §3.3 is 
to note that they scale - differently than in the "maximal merging" model, see Carlberg 
(1996) - like r" 1 = r" 1 (M^/M) 1 / 3 (1 + z) 2 - 5 . This is because in the rate r" 1 = n g EV 
the scaling of the background density (1 + z) 3 is partly offset by the LSS contrast 
C{z) oc (1 + z) -0 - 5 ; with mass, on the other hand, the scaling results from n g oc M _1 
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and E oc M 2 / 3 . Toomre (1977) has fixed the observed interaction probability t /r at 
a few percent for interactions of bright local galaxies; we find r" 1 = C n g (t ) E V a ~ 
1/500 Gyr -1 . At z « 1 our probability £/t for galaxies of 5 x 10 9 Mq exceeds 30 
%. The bright career of this dwarf population activated by interactions is effectively 
started at z 2 when the large scale structures begin to loom out; they reshuffle the 
conditions set by the canonical hierarchical clustering including the initial correlations of 
linear perturbations, sustain contrasts larger than unity and produce non-linear relative 
velocities confined to one or two dimensions. This career terminates, due to the Hubble 
expansion, with a progressive fading for z ^ 0.5 and with the number of sub-L* units 
reduced by ~ 10 _1 . 

The LSS features used in our computation - filamentary (D = 1) or sheet-like 
(D = 2 as considered above) structures produce similar results - come down to three 
key points: formation era around Zi n ~ 2, local density contrast ~ 5, and slow growth 
of the contrast from Z{ n to the present. These features, which appear in the relatively 
local observations and also in all N-body experiments, are also incorporated in analytic 
interpretations different from those of Zel'dovich's (see Bond, Kofman, & Pogosyan 
1996). We stress that all such structures are ultimately traced back to long-distance 
correlations, not considered in the homogeneous renditions of the hierarchical cluster- 
ing scenario, but actually intrinsic, and effective in bringing to interact galaxies that 
otherwise would not have met. 

The detailed value of the initial contrast C is the only "free" parameter at present. 
In fact, its magnitude is constrained by its very definition to be of the order of a few, as 
we have used. As to finer adjustments, its local distribution is currently being measured 
(Ramella, private communication); the analytic theory of its evolution is being refined, 
beyond the approximation used in §3.3 (see Sahni & Coles 1995); observations are 
reaching beyond z«l (see Willinger et al. 1996). 

Our results are conserved when the shape of the initial mass function is changed 
from the detailed Press & Schechter form but is still steep as predicated by all hierar- 
chical clustering scenarios (cf. White 1994). 

The main trends are also robust with respect to variations of the M/L ratio from 
the value M/L oc M" 1 / 3 used in §§4 and 5, to include M/L oc const, and M/L oc M" 2 / 3 
as shown in Figure 5. The different M/L ratios tilt, but moderately, the shape of the 
distant LFs; correspondingly, the tail of N(z) for z > 1 rises to the values indicated by 
the data by Cowie et al. (1996) if the dependence of M/L goes toward M -2 / 3 . 

We also computed the counts in the K-band. These, as expected, increase consid- 
erably less than in the blue due to the long evolutionary times of the red stars and to 
the stronger /c-correction. 

One recurrent objection to the aggregation scenario concerns the correlation func- 
tion, as discussed by Efstathiou (1995). This is a delicate issue, as is shown in Figure 
6. At the statistical level of the correlation function, aggregations with cross section 
less than linear in mass (like E oc M 2 / 3 for nearly geometrical cross sections that ap- 
ply here, see §2) do not cause the slope of the initial correlation function to change. 
A steepening instead takes place for the stronger, gravitationally focussed interactions 
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with non-linear cross section E oc M 4 / 3 . The whole issue may peter out, however, with 
the growing evidence for clustered redshifts at z 0.8 — 1 as observed and discussed by 
Le Fevre et al. (1994) and by Koo et al. (1996). 

The other recurrent objection to aggregations is that the spiral disks may be dam- 
aged by merging, as discussed by Ostriker (1990). In our framework, however, at the 
level of M* the mass growth is due to minor episodes of merging with the numerous 
but much smaller (and often more diffuse) galaxies. On the other hand, the overall 
change of the characteristic mass seen in Figure 4a is mainly due to the changing shape 
of the mass function, i.e., to its flattening corresponding to the disappearance of small 
galaxies. This results in a mild statistical change in time of the second moment of the 
distribution, which does not imply a corresponding change in the masses of individual 
galaxies. For similar reasons, the population of ellipticals does not evolve appreciably 
by aggregations in the redshift range considered here. 

To conclude, excess counts at B ^ 24 may result from the combination of several 
effects: local shape and normalization of the luminosity function, long line-of-sight as- 
sociated with fi G < 1, passive evolution, and interactions as stressed here. The shape of 
N(z) is more telling, and at relatively bright magnitudes mg ^ 24 requires some evo- 
lution, yet it rules out the pure luminosity kind in a flat universe. The resolved N(L, z) 
with its evolving shape restricts the issue to the following alternatives: luminosity- 
dependent luminosity evolution or a combination of (differential) density plus luminos- 
ity evolution, as naturally provided by aggregations. The steepening associated with 
the latter drives a specific differential shift of N(z), characterized by the peak marching 
on toward higher z at fainter magnitudes (see Figure Ad); it also sustains the rise of 
the counts at very faint magnitudes (see Figure 4c) even beyond bj = 27 measured by 
Metcalfe et al. (1995), probing in an integrated fashion the driver's action at redshifts 
not yet reached by spectroscopy. 

The above predictions stem from simple dynamics in which interactions produce 
aggregations and starbursts together. On the other hand, Cavaliere & Menci (1993) 
and Moore et al. (1996) point out that in clusters higher densities but larger relative 
velocities cause more frequent but weaker encounters ("harassment") resulting mainly 
in pure starbursting and in blueing of the Butcher-Oemler (1984) type. In fact, such 
flybys and the true aggregations constitute two facets of a single process, differing as for 
the detailed form of the cross section (see equation 2.3) and for the ratio s(z)/t~ 1 (z) 
of the star formation to the interaction rates. 

These two facets stem from the same dynamics, and concur to stress the importance 
of interactions as drivers of galaxy evolution. Actually, the intermediate regime should 
also be observable in clusters, in the form of a local LF frozen with a steeper faint end 
or a more pronounced upturn, as shown by Figure 3e. 

We have benefited from many informative discussions with G. De Zotti, M. D'Onofrio, 
E. Giallongo, M. Ramella, R. Windhorst and G. Zamorani. We thank the referee for 
his helpful comments which improved our presentation. This work was supported by 
partial grants from MURST and ASI. 
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FIGURE CAPTIONS 



Figure 1: Results for galaxies with M/L oc M x / 3 interacting in the homogeneous 
"field" with = 1 and A = 0. 

Panel a): The solid lines show the comoving mass function at 10 equally spaced 
time intervals, from the initial time ti n (dashed line) to the present. 
Panel b): The corresponding luminosity functions in the B band for three redshift 
bins: z = — 0.2, 0.2 — 0.5 and z > 0.5; the data are taken from Ellis et al. 1996. 
Panel c): The computed galaxy counts (solid line), compared with data from various 
authors: Maddox et al. (1990) solid circles; Jones et al. (1991), open triangles; 
Tyson (1988), open circles; Lilly, Cowie & Gardner (1991), open stars; Metcalfe et 
al. (1995), solid triangles. 

Panel d): The computed redshift distribution in the magnitude range 22.5 < uib < 
24 is compared with the data from Glazebrook et al. (1995). 

Figure 2: Same as Figure 1 but for fi G = 0.2 and A D = 0.8. 

Figure 3: Same as Figure 1 but for galaxies within virialized structures. 

Figure 3e: The effect of varying the position of the cutoff for the efficiency e(V/v m ) in 
the cross section (2.3). In this figure e = is taken for V > 2.5 v m rather than 
V > 3 v m as adopted throughout the paper. 

Figure 4: Same as Figure 1 but for galaxies within large scale structures with D = 2. 
In panel b) we show also our prediction for the luminosity function at z = 1 (dot- 
dashed line). In panel c) we show also our predictions for 24 < tub < 26. 

Figure 5: The evolution of the luminosity function of interacting galaxies within large 
scale structures, computed for M/L =constant (left panel) and for M/L oc M -2 / 3 
(right panel). 

Figure 6: The time evolution of the slope 7 of the two-point correlation function for 
aggregating galaxies, as computed from a typical Montecarlo simulation (see Menci, 
Colafrancesco, & Biferale 1993 for details); 7j n is the initial value. The dotted line 
shows the result for a cross section scaling as S ~ M 4 / 3 , and the solid line that for 
E ~ M 2 / 3 . 
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Fig. 6 
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Table 1: The parameters for the aggregation rate given in equation (2.7) and used in the computation are 
tabulated for the different galaxy environments discussed in the text; Zi„ is the starting redshift, Vi n is the 
initial galaxy relative velocity, C is the initial effective density contrast in each environment, and f(z) is the 
evolutionary factor of the aggregation rate. D defines the dimensionality of the LSS where galaxies reside (D = 1 
for filaments, D = 2 for sheet-like structures). 



